P-wave pairing and ferromagnetism in the metal-insulator transition in two 
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Based on recent experimental evidence for a spin polarized ground state in the insulating phase 
of the two-dimensional electron system, we propose that ferromagnetic spin fluctuations lead to 
an attractive interaction in the triplet channel and cause p-wave pairing in the conducting phase. 
We use the Landau Fermi liquid phenomenology to explain how the enhanced spin susceptibility 
near the critical density yields an attractive potential, in a similar mechanism to superfluidity 
in ^He. As the density is decreased, the p-wave order parameter undergoes a transition from a 
unitary to a nonunitary state, in which it coexists with ferromagnetism for a range of densities. 
As the density is further reduced, the pairing amplitude vanishes and the system is described by a 
ferromagnetic insulator. Thus, we find two quantum critical points as a function of density associated 
with the polarization of the paired state and ferromagnetism. We explain the magnetotransport 
measurements in parallel and perpendicular magnetic fields and propose a shot noise experiment to 
measure the pair charge. 



PACS: 71.10.-w, 71.30.-|-h, 73.23.-b 



I. INTRODUCTION 

It has become self-evident that the spin properties of 
interacting electrons or holes in two-dimensions (2D) play 
a central role in the transport properties of these systems, 
in particular in the possible metal-to-insulator transition 
(MIT) observed in a number of different material sys- 
tems, such-as Si-MOSFETs, n- and p-doped GaAs, AlAs, 
and SiGclil Recent data in Si-MOSFETs by Shashkin et 
al.a and Vitkalov et alB on the saturation of the conduc- 
tance as a function of the magnetic field parallel to the 
2D plane (see also Ref. ^) , combined with previous anal- 
ysis p£-Shubnikov-de Haas oscillations in tilted magnetic 
field,Eli3 suggest that the insulating state is spin polarized. 

These experimental observations revive an unresolved 
theoretical problem on the possible phases of electronic 
systems in 2D as a function of the interaction strength 
(or, alternatively, the densitjd even in the idealized clean 
systems. Bergman and RiceQ raised the possibility that 
as the density is decreased, there is a transition from 
a paramagnetic Fermi liquid into a ferromagnetic Fermi 
liquid state. |Quantum Monte Carlo studies by Tanatar 
and Ceperleyu have considered three different electronic 
states, a paramagnetic liquid, a ferromagnetic liquid, and 
a Wigner crystal, and found a transition from the para- 
magnetic liquid to the Wigner crystal at r^ « 37. How- 
ever, the energies of these three phases become rather 
close for a range of r^. Thus, it is not unreasonable that 
either improved energy estimations, or disorder effects, 
may bring the energy of the ferromagnetic state to lower 
values, so that it may exist for a window of densities be- 
tween the paramagnetic liquid and the Wigner crystal. 
Indeed, perturbative Renormalization Group (RG) cal- 
culations for disordered and interacting electrons in 2D 



by FinkelsteinEl have pointed out a runaway flow in the 
triplet channel even in the limit of low densities (see also 
Ref. no) . This runaway has been recently interpreted as a 
tendency towards ferromagnetism in the diffusive (metal- 
lic) regime,OE3 suggesting that disorder may also help 
trigger a full spin polarization at low densities. Recent 
numerical studies have also indicated a strong tesidency 
towards spin polarization in the localized regimeil3 
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FIG. 1. Variation of the relevant order parameters as a 
function of electron density (mean-field theory), m is the 
magnetization, \l/ is the p-wave pairing amplitude, and 5* de- 
notes the pair spin projected along m. Spontaneous mag- 
netization sets in at the critical density no, while at ric the 
systems becomes a disordered ferromagnetic insulator, ripoi 
is the density where the pairs become fully spin polarized. 



Provided the recent experimental evidence from trans- 
port measurements (not susceptibility measurements) for 
the ferromagnetic state of dilute 2D electrons, we in- 
vestigate the effect of enhanced spin fluctuations on the 
paramagnetic side of the transition. Wc assume that the 
paramagnetic state can be described, for densities above 
a critical density no and for finite temperatures below 
the Fermi energy, by the Fermi liquid phenomenology. 
As one approaches the critical density, one of the Lan- 
dau Fermi liquid parameters, Fq, which renormalizes the 
spin susceptibility, crosses the miiiMnum bound for the 
Pomeranchuk's stability condition.llil We argue that the 
proximity to the ferromagnetic instability leads to an at- 
tractive interaction for a range of densities preceding the 
ferromagnetic transition. The attraction, in the triplet 
channel, leads to p-wave pairing. Once the system en- 
ters the ferromagnetic state, the paramagnon exchange 
mechanism for the attraction rapidly decreases, and these 
two phases compete to the point where the p-wave state 
ceases to exist at a density ric- Therefore, there are two 
quantum critical points in the phase diagram at densities 
ric and uq (likely very close) , and there is an intermediate 
region in densities where the two phases coexist, but in 
which the p-wave pairing is in a nonunitary state. At the 
mean field level, the phase diagram is shown in Fig. n^. 

The possibility of unconventional pairing in an elec- 
tronic system has been recently considered in the con- 
text of layered Sr2Ru04. In the Sr„+iRu„03„+i series, 
where n determines the number of Ru02 planes in the 
unit cell, ferromagnetic states are observed for n > 3. It 
is believed that the proximity to ferromagnetism plays 
an important role for superconductivity in the Ruthen- 
ates, and p-wave spmajetry was proposed on the basis of 
similarities to '^He.Ej'EEl Although a direct transition be- 
tween ferromagnetic and superconducting states has not 
been experimentally observed, the Ruthenates show sim- 
ilarities to the problem studied in this work. We should 
point out, however, that in Ruthenates there is true 3D 
long range order, while we discuss systems that are truly 
two-dimensional and therefore subject to strong fluctua- 
tions. 

We also draw a strong analogy between the 2D elec- 
tronic-states and the 3D '^He systems in a number of 
ways.cZI For one, we argue that the Fermi liquid phe- 
nomenology should not be dismissed in describing the 
2D interacting electronic system in the metallic side of 
the transition at finite temperatures. One of the usual 
concerns that is raised against the Fermi liquid state in 
the 2D problem at low densities is that the ratio between 
the Coulomb and kinetic energies is about a factor of ten. 
This logic can be misleading, a) since the contributions 
from exchange and correlations reduce this ratio, and b) 
Fermi liquid parameters, which measure the pirength of 
the interaction relative to the kinetic energy,Ej are typ- 
ically large compared to unity even for 3D ■^He. For 
example, the Landau parameters for '^He at high pres- 
sure (27 bar) that renormalize the compressibility, mag- 
netic susceptibility, effective mass, and spin precession 



rate are F^ 



68.17, F;} = -0.76, Ff = 12.79, and 
Ff — — I.OO.E2 Even though most Landau parameters 
are very large when compared to unit, the normal phase 
of '^ He is very well described by Landau's Fermi liquid 
theory. Indeed, if any, the real question is why Landau's 
phenomenology works so well, way beyond the perturba- 
tive regime where RG arguments for fermionslij justify 
the stability of the Fermi liquid. 

We show that we can consistently interpret the re- 
cent data Qn-jthe 2D MIT transition close to the criti- 
cal densitymj as due to the enhancement of the Landau 
parameter Fq. The paramagnon exchange mechanism 
can be responsible for an attraction in the triplet chan- 
nel in the 2D electron problem, in complete analogy with 
the problem of 3D "^He. The major difference between 
these two systems is the dimensionality: in 2D true su- 
perconductivity or superfluidity is only possible at zero 
temperature. Strong fluctuations in 2D do not allow for 
long range order. For a singlet paired state in 2D with 
U(l) symmetry algebraic order can be established below 
the Kosterlitz-Thouless (KT) temperaturcEj However, 
for a triplet paired state where the order parameter is 
a complex vector and the symmetry group is nonabeliap- 
superconductivity can only be established at T = O.EJ 
We thus propose that the insulating state is ferromag- 
netic, while the metallic state corresponds to a paired 
p-wave state. The possibility of singlet superconductiv- 
ity in the observed conducting phase Q£-the 2D electron 
systems waSj^iggested by Phillips et alE2 and Belitz and 
Kirkpatrick;E3 here we present a possible mechanism for 
pairing (without finite temperature long-range order) in 
the triplet channelo 

We would like to stress that we do not address in this 
paper the reason why the measured conductance in Si- 
MOSFETs seems to saturate in the triplet paired state 
when T ^ 0. In the case of singlet paired superconduct- 
ing 2D materials (thin films) and 2D Josephson junc- 
tion arrays, where a finite KT transition should be ob- 
served together with the vanishing of the resistpity, a 
saturation of the conductivity is also observed. E3 The 
source and precise mechanism for these dissipative effects 
are preae|Li±ly unknown, despite some recent theoretical 
efforts.OEZl 

There is a natural question regarding the possibility 
of a paired state in the 2D electronic systems where the 
MIT is observed: Why can one have pairing if the con- 
ductance is of order e^//i near the transition? Naively, 
if one uses intuition from noninteracting electrons, then 
the bare value of the conductance is Go = {2e'^/h){kF£); 
weak localization corrections, perturbative in {kp£)^^, 
are added to this bare value. This would imply that 
kpl ~ 1 at the transition, and therefore disorder is too 
strong and pair breaking. However, the value of £ that 
one reads from this naive argument is not a measure of 
disorder alone. The energy scale of the interactions is 
larger than the Fermi energy of the 2D electron systems 
near the transition, thus the scattering of electrons even 
if the disorder is weak should be large and dominant. 



Therefore, a dimensionless conductance of order unity 
does not necessarily imply strong disorder. For example, 
it is known that near 2D superconducting-insulator (SC- 
I) transitions the dimensionless coiid|Uctance is of order 
one even in the absence of disorder .E3 We claim that the 
same happens in the context of the MIT transition dis- 
cussed here. We also give arguments showing, based on 
Fermi liquid theory, that near the ferromagnetic transi- 
tion the paramagnetic scattering can provide for conduc- 
tances of order unity. 

The paper is organized as follows: In Section H we dis- 
cuss the Landau's Fermi liquid phenomenology applied to 
the paramagnetic phase of the 2D electronic system, and 
the Pomeranchuk instability leading to ferromagnetism 
at low densities. We argue that the enhancement in spin 
fluctuations due to the proximity to the ferromagnetic 
state causes an attractive intera ctio n in the triplet chan- 
nel and p-wave pairing. Section III contains a discussion 
of a two-component model for the MIT transition, where 
we argue that, due to electron-electron interactions in 
the presence of disorder, the dimensionless conductance 
due to paramagnetic scattering can become of order unit 
close to the transition. The mean-field phase diagram of 
the problem is established in Section IV using a generic 



Ginzburg-Landau free energy for a p-wave paired state 
coupled to a ferromagnetic order parameter. In Section 
M we compare our results to the available experimental 
data for parallel and perpendicular magnetic fields for 
the 2D electron gas in a Si-MOSFET and propose new 
shot noise experiments that can test our theory. Section 
VI contains our conclusions. 



II. THE P-WAVE PAIRED METALLIC PHASE 

Our starting point is a Landau Fermi liquid theory for 
the metallic phase of the 2D electronic system. We con- 
sider first the effects of interactions, and then those of 
disorder. Let us briefly review and then apply the Lan- 
dau phenomenology to the experimental observations on 
Si-MOSFETs. Consider an isotropic Ferm i liquid with 
planar density n, Fermi momentum kp = ^Jli^nj g^^ and 
Fermi energy Ep — h kp/2m* — irfi n/g^m* , where m* 
is the effective mass and g^ accounts for the valley de- 
generacy. The ground state of the problem is described 
in terms of quasiparticles that fill up a Fermi sea up to 
the Fermi energy. The change in the energy of a Fermi 
liquid due to changes in the quasiparticle^iflxge density, 
6n{k), and spin density 5a{k) is given bylUil3 

5E= f (fkei^Sn{k)+ f tfktfk' f{k,k') Sn{k) Sn{k') 



(fk cPk' r{k, k') 5a{k) ■ Sa{k'), 



(1) 



where e^ is the bare dispersion and /*(/c, k') and /"(A:, k') 
are the symmetric and anti-symmetric Landau parame- 



ters, respectively, 
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In 2D these parameters can be ex- 
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where n gives the angular momentum in the plane, and 
9^^:, is the angle between k and k' . It is useful to de- 
fine dimensionless parameters f"'" = Af(0)/*'", where 
N{0) = g^m* /ttTl is the 2D density of states at the Fermi 
energy. 

The stability of the Fermi liquid state (or the Fermi 
surface) is given— in. Landau's theory, by the Pomer- 
anchuk criterionJiJH which in 2D can be written as 



>-l 



(3) 



for all values of n. Since all the physical quantities in 
Landau's theory can be written in terms of the Landau 
parameters, a violation of the Pomeranchuk criterion im- 
plies an instability of a physical observable. The com- 
pressibility, for instance, is given by 



1 + Fn^ 



(4) 



andf-an instability to phase separation implies that F^ < 
— 1.E2I In the same theory the effective mass is given by 
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(5) 



vanishing when F^ = —1 (here rrii, is the carrier band 
mass). The magnetic susceptibility can be written as 
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where gg ~ 2 is the bare (band) Lande g- factor. Thus, for 
Fq w — 1 the magnetic susceptibility diverges, indicating 
an instability towards a magnetically ordered phase. 

As one of the Landau parameters approaches the crit- 
ical value given by Pomeranchuk's criterion, there is a 
strong enhancement of the interactions in the Fermi liq- 
uid. Consider, for instance, the case of density-density 
interactions that are determined by Fq. It can be shown 
that the induced density-density interaction in a Fermi 
liquid (in pUae static limit) is given by the usual RPA 
expressiorH 
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Thus, close to the Pomeranchuk's instability Fq 
the interaction is very large and attractive, leading to 
phase separation. On the other hand, the induced spin- 
spin interactions in the same system is given by 
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(8) 



where a is the electron spin. Therefore, when F^ w — 1, 
that is, close to the magnetic instability, this interaction 
is also large and attractive, leading to pairing in a spin 
triplet channel. In fact, we can estimate the size of the 
pairing amplitude using the weak coupling BCS expres- 
sion 



Ef e 



(9) 



where the Fermi energy Ep works as a cutoflLin the prob- 
lem because is the only energy scale present.E2l Obviously, 
we have to consider Eq. (^) carefully since, as the system 
approaches the instability, the attraction is very strong 
and the weak coupling expression breaks down; in this, 
case, one should use a strong coupling approximation.Eil 
Thus, the expression in Eq. (O) can still be used when the 
attraction is weak \F^\ <^ 1 and is only a crude estimate 
when Fg Ki-l. 

In general, one expects the Landau parameters to be 
dependent on the electronic density n. Let us consider 
the situation of a Fermi liquid close to a magnetic insta- 
bility that happens at n = riQ. The Landau parameter, 
Fq{5)^ can be expanded close to the transition as 



where a > is a constant and 
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(10) 



(11) 



measures the distance from the quantum critical point. 
In Eq. (|l^) we disregard higher order terms in the den- 
sity variations around the critical point. Observe that in 
this case the magnetic susceptibility can be written from 
Eq. (PI) as 
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(12) 



showing that the susceptibility diverges linearly with the 
distance from the critical point. Notice, from Eq. (g), 
that the weak coupling expression for the pairing ampli- 
tude close to the transition is given by 



|Ap(5)| « Epe 
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(13) 



This result indicates that the pairing amplitude is of or- 
der of the Fermi energy in the system. In the case of the 
2D Si-MOSFETs, Ep is usually of order of 5 K because 
of the low electronic density. The critical temperature, 
Tc, however, remains zero because of the dimensionality 
of the system, otherwise this problem would be a case 
of high temperature superconductivity. Notice that the 
number of pairs is Up = n/2 at zero temperature. 

Away from the transition we can classify the behavior 
of the system depending on the full density dependence 
of Fq. At large enough densities, n > ha, Fq should be- 
come positive because of the screening of the electronic 



interactions by backflow effects. In this case, without the 
help of pairing, the localization effects should dominate, 
and the system should become an Anderson insulator. 
So it is only in this high density (n > ha) regime that 
one could possibly attempt to apply ideas established 
through the scaling theory of localizationE3 for noninter- 
acting electrons. For tiq < n < ua we have — 1 < Fq < 
and pairing is effective in delocalizing the electrons, lead- 
ing to the metallic state observed experimentally. In this 
regime, naive intuition based on noninteracting electrons 
should not apply. Finally, for n < uq, the system be- 
comes a ferromagnet. We will return to this regime later 
when we consider a Landau-Ginzburg theory describing 
the system. We show that there is a coexistence region 
with both p-wave pairing and ferromagnetism for den- 
sities Uc < n < no, and for n < Uc the physics of the 
problem is the one of a ferromagnetic insulator. 



A. The Effect of Disorder in the Paired State 

We have argued that the proximity to the ferromag- 
netic instability induces pairing in the Fermi liquid state. 
The effects of strong interactions are built in in this pic- 
ture, but we still need to discuss the effect of disorder in 



this paired state (which we argue in Section IV should be 
in fact a fully gapped px ± ipy state) . The question to be 
addressed is: will the paired state survive disorder? As 
it is known for the case of strongly coupled superconduc- 
tors, which is the case here, when the electron mean free 
path, £, becomes of the acjder of the coherence length, 
^sc, pairing is suppressed.EJ Since the pairing amplitude 
is essentially of order of Ep [see Eq. (Oh], the coherence 
length is of the order of the Fermi wavelength. 
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and of the size of the intcrparticle spacing. Therefore, 
pairing should survive as long as i* 3> ■^sc, or equivalently, 
kp£ 3> 1. The main question here is whether this condi- 
tion is satisfied in the heterostrutures where the MIT is 
observed. We argue below in favor of this case. 

The common belief is that close to the MIT transition 
kpi « 1 because the dimensionless conductance at the 
transition is of order one. This indeed would be the case 
if the value of the conductance would be completely fixed 
by the amount of disorder in the sample. However, as we 
have shown, there are strong electron-electron interac- 
tions close to the ferromagnetic transition (indeed inter- 
actions are the reason for the instability, since Fq — > — 1). 
These interactions in the presence of localized electronic 
spins, as we are going to show in Section HI, can provide 
a large contribution to the resistance. In other words, the 
value of the critical resistivity, pc, can be of order h/e'^ 
from interaction effects even if disorder is small. The 
concept of using {kp£)~^ as a measure disorder is only 
good if one could determine what £ is independently of a 



conductivity measurement, because otherwise strong in- 
teraction effects will blend in and make such discussion 
useless. Therefore, in strongly correlated systems one 
should be careful in extending arguments that are only 
valid for noninteracting systems. 

/^From a theoretical perspective, it is known that in- 
teractions alone can lead to a univeriial conductivity of 
order ctq = 4e^/ft.. It has been showncj that, for the 2D 
superconductor-to-insulator transition in the clean Bose- 
Hubbard model (belonging to the 3D-XY model univer- 
sality class), the conductivity is given by cr « 0.285cro. It 
was argued in Ref. ^ that this transition is in the same 
universality class of a superfluid-insulator transition for 
bosons moving in a random potential. In 1 + 1 dimen- 
sions, it is known that the universality class of a super- 
conductor to insulator transition for a model of fermions 
with attractive interactions is the same as the insulator 
to superfluiiLtransition in a model of repulsively interact- 
ing bosons.O Since 1-1-1 is-the lower critical dimension 
for this type of transition,Ea we expect this type of re- 
sult to hold in our case as well. Indeed, there is strong 
experimental evidence for these results in the supercon- 
ducting to insulatpi; transition in amorphous supercon- 
ducting thin films E3 



III. THE TWO-COMPONENT MODEL 

Based on experimental evidence that the insulating 
phase in the 2D electron system is a ferromagnetic state, 
we have argued in Section |l^ that the proximity to the 
spin polarized state induces p-wave pairing on the metal- 
lic side of the MIT. In this Section we study the transition 
between the p-wave paired state and the ferromagnetic 
state using a two-component model. Let us split the total 
electron density at T = into a localized ferromagnetic 
component, n^, and an itinerant, paired component, ?!„: 



n — TiL + 2np 



(15) 



The dependence of both n^ and np as a function of n 
is represented schematically in Fig. 0. For densities 
no < n < ua, all electrons are paired {np — n/2) at 
T = 0, while for very low densities, all electrons are in 
the insulating ferromagnetic state {n = ul). There is a 
region where the two components coexist: as the system 
starts to spin polarize below ng , the density of pairs starts 
to decrease, vanishing at some critical density ric < no- 
At very low densities (n ^ Uc) the Wigner (rrystal or 
Bragg glass phase becomes the ground state.HE3 In this 
work, however, we do not attempt to study such phase, 
since it requires unscreened long range interactions that 
are not included in our formalism. The study of the 
Wigner crystal would then need a different starting point, 
far from the MIT. For this reason, the Wigner crystal is 
absent in the phase diagram of Fig. ||. 

The appearance of a coexistence region is natural in the 
way we divide the two densities. It also follows from the 



competitive nature of the two phases, as the mechanism 
that leads to pairing - paramagnon exchange - is strongly 
suppressed as the system spin polarizes. The coexis- 
tence phase is also found within the Landau-Ginzburg 
theory of Section IV, as a direct consequence of having 
two competing orders. The existence of metallic regions 
surrounding localized puddles of charge has been con- 
firmed experimentally through ingenious local compress- 
ibility measurements in GaAs/AlGaAs samples.EZl These 
experiments reveal that the fragmented localized regions 
increase in number as the density is lowered towards the 
MIT. The appearance of these fragments in the metal- 
lic phase (n > nc, no) is consistent with the idea that 
some incipient magnetism is always present in diffusive 
systemS[-even when the interaction in the triplet channel 
is weak.EHI In regions where strong magnetization occurs 
due to fluctuations in the wave functions, the param- 
agnon mechanism is absent, thus suppressing pairing and 
leading to localization. This is not taken into account in 
our mean-field (homogeneous) treatment of the problem; 
had it been considered, it would likely prolongate the n^ 
curve towards densities larger than no (not shown in Fig. 




FIG. 2. Densities of the localized, ul, and itinerant, 2np, 
components as a function of total density near the critical 
region. 

Next, we use this two-component picture to provide an 
argument on why the critical conductivity near the MIT 
is of order e^/h. Consider the situation at finite tem- 
perature within the coexistence window Uc < n < uq. 
While at zero temperature all the delocalized electrons 
(with density 2np) participate in the pairing, at finite 
temperatures only a small fraction is really paired. The 
unpaired state can be described by the Fermi liquid the- 
ory of Section O. The localized ferromagnetic component 
is a source of spin scattering for the itinerant Fermi liq- 
uid. The coupling 



SH^J d\ SL{r) ■ Sitir) 



(16) 



follows from splitting the electrons into two components 
in the effective Hamiltonian for the interacting electron 



liquid. Here ^'^(r) and Sit{r) are the electron spin oper- 
ators for the localized and itinerant components of the 
electron system, respectively. The exchange coupling 
J = iV(0) Fq can be obtained from the expression for 
the energy in the Landau Fermi liquid phenomenology. 

This interaction leads to a high temperature scattering 
time for the itinerant fluid given by 
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for scattering off the spin S = 1/2 ferromagnetic compo- 
nent. This gives a resistance 
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-^{F,rSiS + l). (18) 



For a narrow coexistence region, the Landau parameter 
Fq should be close to the Pomeranchuk critical value of 
— 1; hence, the conductance is solely determined by the 
ratio X — nL/n: 
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Notice that this result implies that there can be a large 
variation of the resistivity in the critical region of the 
phase diagram. The scale for resistance is the prefactor 
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(20) 



and if near the MIT transition region the fractions of the 
two components are close, the ratio x/{l — x) should be 
of order unit. Hence, the high temperature resistance 
near the "separatrix" line should be of order pg, which 
gives a conductivity of order a ~ 0.3 ctq- This simple 
argument highlights the importance of considering the 
electron-electron interactions when calculating the resis- 
tivity of the system near the ferromagnetic instability 
(^0 = -!)• 



IV. LANDAU FREE ENERGY - MEAN-FIELD 
PHASE DIAGRAM 

In this section we discuss the transition between ferro- 
magnetism and p-wave superconductivity. In the context 
of Ruthenates it has been shown that the ferromagnetic 
to superconductor tpjpaition can be described in terms 
of an SO(IO) modeOo Here we do not take the high 
symmetry approach; instead, we simply write a Landau- 
Ginzburg free energy that combines both p-wave super- 
conductivity and ferirOmagnetism (we assume the system 
to be homogeneous) :E3 



F[m, ip, d, d* 
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— [zd(0)xd*(0)].m, (21) 

where m is the ferromagnetic order parameter and ip and 
d{9) are amplitude and vector parts used to describe the 
p-wave pairing order parameter at the Fermi surface: 



*a/3(6') = iipy^}crkcr2]af3dki0) 



(22) 
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with the 3D vector d{9) obeying the normalization con- 
dition 
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(23) 



Here ak with k = 1,2,3 are Pauli matrices. We assume 
that the coefficients S and 7 appearing in the free en- 
ergy are nearly independent of electron density n, while 
a = a{n — ua) and b = f3(n — ns)- All a, (3, 7, and 5 are 
positive. 

Clearly not all fourth order terms allowed by symmetry 
have been taken into account in Eq. (|2|). The inclusion 
of all terms would render the analysis extremely difficult 
even at the mean-field level. Thus, the choice manifest in 
Eq. (^ should be considered as the simplest one that 
reproduces the phases discussed in previous sections. 

It is worth noticing that the expectation value of the 
Cooper pair total spin operator at a point of the Fermi 
surface is given by 



(S) ^iiP^diO) xd*{0). 



(24) 



Thus when d{9) is a real vector, apart from a overall 
phase factor, (S) = 0. In this case, named unitary, one 
can show that d{9) defines a direction along which the 
spin operator S has eigenvalue zero. 

One could in principle minimize the free energy with 
respect to all ten real parameters, namely, Tp, rrii, di, 
d*, with i = 1,2,3, using a Lagrange multiplier to en- 
force the normalization condition. Instead, we restrict 
the form of the vector d{9X-ijO certain classes, following 
the treatment used for ^He.l23 For instance, the analogous 
to the B phase or Balian-Werthamer (BW) phase of '^He 
would correspond to the isotropic (nodeless) choice 



d{9) = ei cos 9 + 62 sin9, 



(25) 



which is clearly unitary. Another possible unitary choice 
(also nodeless), resembling the A phase or Anderson- 
Brinkman-Morel (ABM) phase of "^He is 



d{9) — 63 (cos 9 ±i sin t 



(26) 



(The latter case seems 
ing phase of Sr2Ru04 



f relevant to the superconduct- 
5l|i6|'j Here, we use neither choice. 



but adopt instead a parameterization that allows for a 
nonunitary d{9), thus (S) 7^ 0. Recall that there exists 
a ferromagnetic coupling between the local magnetiza- 
tion and the Cooper pair total spin (the last term in the 
free energy). A unitary order parameter would make this 
coupling to vanish identically. Our choice will be 



(*a./3) = 



(35) 



d(0) 



V2 



(wei +w* 62) , 



(27) 



where z ; 



and |wp — 1. One may ask how this comes 



about. The reasoning is simple. First, we look for a mini- 
mum energy configuration and therefore fix it to be node- 
less. Second, for a 2D p-wave paired state, we must have 
di{h) = DiaUa, with i = 1,2,3, and n = (cos0,sin0). 
Thus, we can write that A{9) = z vi -I- 2;* V2, with Vi and 
V2 two complex three-dimensional vectors that do not 
depend on 6, but must obey the relation jvip -|- |v2p = 1 
due to the normalization condition on d{0). It is straight- 
forward to check that the BW and ABM states corre- 
spond to the choices Vi_2 — (ei =F «e2)/-\/2 and Vi = 63, 
V2 = 0, respectively. On the other hand, our choice of pa- 
rameterization can be obtained by setting, for instance, 
vi = [w ei + w* €.2) 1^/2 and V2 = 0. The choice is not 
unique and has to be considered as an Ansatz. 

We can now proceed with the minimization of the free 
energy. With our choice for d(0), wc find that 



|d(0)|^ 



1 



identically 


and 
id 


'9) X 


d*{e) = - 


An (pes, 
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gj0/2 


. As 


a result. 
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- aip' 
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sin^ ip 
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m 



(28) 



(29) 



-\-5 1pm + jtp rris sin (j). 



(30) 



(Notice that we do not need to use a Lagrange multiplier 
to enforce the normalization condition anymore.) The 
minimization equations read 



= a + -0^ [1 -I- sin^ 0] + i5 ^n^ + 7 "t-3 sir 
= mi,2 {b + rn^ + 25 ip'^) , 
= ms (5 + m^ + 25 ^'^) +71/'^ sin 0, 
= V' cos (p [tp sin 0-1-7 ms\ . 



(31) 
(32) 

(33) 
(34) 



After Eq. (|32|), we can set mi, 2 = without loss of 
generality (thus m — mj, hereafter). Equation (p4|), how- 
ever, is satisfied only when 'p'^ — 0, or cos(t> — 0, or 
m = —'pP' 1^ sin0. In our analysis of the possible solu- 
tions to these equations we will be looking for the fol- 
lowing sequence of phases as the electron density is low- 
ered: metal/p-wave pairing/p-wave pairing -I- ferromag- 
nct/fcrromagnct. 

We should stress that because of the dimensionality 
and the symmetry of the order parameter we always have 



at any finite temperature, although V' 7^ ^tnd (S) ^ 0. 
This is possible because, on average, (d(ff)) — and 
therefore Eq. (pSJ) follows directly from Eq. (p2|). 



A. P-wave Paired Phase 

Let us assume that ua 3> ub- Thus, starting from a 
paramagnetic phase and lowering the density, we first en- 
counter a second order phase transition to a paired state 
("0 7^ 0) at n — UA- Indeed, setting ?7i = 0, we find from 
Eqs. ^ 



and (34) that sini/) — and 
ip = —a — a{nA — n). 



(36) 



Presumably, in Si-MOSFETs, ua is a high density, out- 
side the range explored in the experiments that probe 
the MIT. Notice that since sin0 = 0, (S) = and this 
p-wave paired phase is unitary. 



B. P-w^ave Paired + Ferromagnetic Phase 
Coexistence 

Since ub ^ «a, we expect ferromagnetism to appear 
at much lower densities (ns should close to the critical 
density of the MIT). In this range, we may simple take 
a ~ —auA- Looking for a solution with ■0^ 7^ 0, m 7^ 0, 
and sin0 7^ simultaneously, we obtain from Eq. (B4) 
that 



m 



(37) 



Now (S) = 7771 63 7^ and consequently the pairing is 
nonunitary. The Cooper pair spin points along the di- 
rection of spontaneous ferromagnetization. Solving the 
other two equations for m and ip'^, we find 



(3{no - n) 



and 



respectively, with 



and 
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(38) 



(39) 



(40) 



(41) 



Notice that we need 5 < l/\/2 for stability. It is simple 
to see that m and ip go continuously to and —anA, 
respectively, as n ^ no- The transition is again of sec- 
ond order and ng represents the critical density where 
the ferromagnetism sets. The pairing order parameter ip 
decreases as one crosses ng. 



C. Full Cooper Pair Spin Polarization 

Once ferromagnetism appears, the Cooper pair total 
spin tends to align itself with the ferromagnetic order pa- 
rameter, m. Full polarization of the pair spin, (S) — l-es, 
occurs when sin0 — —\ and the onset is marked by the 
point where V'^ = ^"m. Combining this relation with 
Eqs. ( p8| ) and ( |39| ) we can find at which electron density 
one reaches full spin polarization by solving the following 
equation for rtpoi 



^pol 




(42) 



It is clear that the order n\ < Up^i < uq is obeyed, con- 
straining full pair spin polarization to occur while ip is 
still finite. For electron densities lower than ripoi, the 
pairing order parameter continues to drop, while the fer- 
romagnetic order parameter grows. It is important to 
remark that although all Cooper pairs are spin polar- 
ized, not all electrons in the system are spin polarized. 
As discussed in Section III , there also exist unpaired elec- 
trons in the whole range n,. < n < hq. In the coexistence 
region ric < n < ripoi, Eqs. ( pl| ) and ( [3^ ) are reduced to 



and 



a + 2^'^ + dm^ --fm^O (43) 



m{b + m^ + 2S'iIj^)-jtP^ ^0. (44) 



D. Ferromagnetism 

As the density is lowered further than ripoi, the p-wave 
pairing parameter drops and the magnetization (con- 
nected to the spin polarization of unpaired electrons) 
grows. At the point where the pairing vanishes, Eq. (E3) 
yields rri^ = —b. The critical density Uc can then be 
found by solving Eq. (H3), namely. 



ub 



an A 



7 



SV^ 



ve 



riB 



(45) 



For n < Uc, we have m = \J —h and the magnetization 
may increase up to its limit value. Obviously, this mean- 
field treatment does not take into account quantum fluc- 
tuations or the interplay between spin interactions and 
the Anderson localization. Whether these effects happen 
above or below Uc will depend on the microscopic details 
which are outside the scope of the Landau- Ginzburg phe- 
nomenology. 

The exact form of the mean-field phase diagram re- 
sulting from the minimization of the free energy depends 
on the values of the coefficients a, /3, 5, and 7. Besides 
the constraint that all coefficients should be positive, it 



is necessary the inequality ub > an a/ 5(3 to hold in or- 
der to Uc be positive (more specifically, ni < Uc < ripoi 
under this condition). This case is illustrate in Fig. [|. 
Moreover, if 7^ « 2a6nA, and both S and a are suffi- 
ciently small, it is possible to have the (experimentally 
accessible) critical densities Uc, ripoi, and uq very close to 
each other. 



V. COMPARISON WITH EXPERIMENTS 

In this section we analyze some of the recent experi- 
ments that probe the magnetic field dependence of the 
conductivity, and connect them to the theory we develop 
in this paper. Let us look separately into the parallel and 
perpendicular field experiments. 



A. Parallel Magnetic Field: Magnetotransport 

Recent experiments have probed the dependence of the 
conductance iuj^i^MOSFETs as a function of an in-plane 
magnetic field. ETcI The experiments show that the con- 
ductance saturates beyond a field iJsat- The saturation 
value is interpreted as the field needed to fully spin po- 
larize the 2D electron system. Further evidence for full 
spin polarization is provided by an exact doubling of the 
period of the Shubnikov-de Haas oscillations when the 
experiment is done in a slightly tilted [magnetic field and 
the in- plane component exceeds iJgat-™ 

In the experiments of Ref. 0, it was observed that the 
magnetic field iJgat (Be in the notation of Ref. 0) needed 
to saturate the magnetoconductance at low temperatures 
scaled linearly with the difference n—ric over a wide range 
of densities. Close to saturation, the scaled magnetocon- 
ductance curves taken at different densities followed the 
same function of the ratio between H/Hsa.t{n). 

In Ref. 0, the magnetoconductance was also found to 
saturate beyond a value of magnetic field iJsat that de- 
pended on the density. Near the MIT, the conductivity 
followed the scaling 



a{H,T,n) -o■sat(?^) = / 



H 



(46) 



Hsat{n,T) ^ 

Again, similarly to Ref. y, i?sat (x n — no for a wide range 
of densities at low temperatures. However, very close to 
the transition, the density and temperature dependent 
saturation field Hgat extracted from the data behaved 
approximately as 



i/,at(n,T) = A{n)^A{nf+T^. 



(47) 



The parameter A(n) is weakly dependent on density, be- 
ing almost constant for a range of densities, and increas- 
ing by about 20% near the critical density no- The pa- 
rameter A(n) is fitted to a form 



A(n) = Ao (n - Tio)'' 



(48) 



with Ao « 2.27 and 7 « 0.6. 

Let us now explain some of these experimental results 
using the theoretical framework we propose in this paper. 
We argued in Section || that, despite the large interac- 
tions, we can still apply the Landau Fermi liquid phe- 
nomenology in order to understand both the instability 
towards a ferromagnetic state and its precursor p-wave 
paired state. One of the Landau Fermi liquid parame- 
ters, Fq, crosses the Pomeranchuk's stability boundary, 
Fq = —1, leading to the ferromagnetic instability of the 
2D electron system. Near the critical density no, we per- 
formed a Landau expansion for the Landau Fermi liquid 
parameter as in Eq. (|o|) . In Ref . Q the degree of polar- 
ization ^ = g* fiBHn/2Ep, where 



harder to obtain. For example, this may be the cause for 
the density dependent prefactor A{n) in Eq. (|47|). 

A general feature of the scaling Eq. (^^ that we can 
explain, however, is why the finite temperature correc- 
tions are quadratic in T for T ^ A and linear in T for 
T S> A. The low-T behavior follows simply from the 
finite-temperature Fermi liquid susceptibility. The high 
temperature behavior, on the other hand, is due to the 
fact that the susceptibility should, in this temperature 
regime, obey a Curie law. 



B. Perpendicular Magnetic Field: Quantum Hall 
Effect 



90 



1 + F- 



(49) 



is the effective Lande 5-factor re 
dau Fermi liquid parameter FqX 



normalized by the Lan- 
4 Thus, 



? 



1 



gof^sHii 



1 + F^ 2Ef 



(50) 



When the field iJ|| is sufficiently strong to fully spin po- 
larize the system, we have ^ = 1. If we neglect nonlinear 
terms in the susceptibility (which in principle could be 
important when the spin polarization is large), the satu- 
ration field in the metallic phase can be estimated from 
the linear response expression of Eq. (pOh , 



ffsat = ^^ (1 + F,'^) 

5oMs 



2Trh^ 



gom^^B 



■aJn-no), (51) 



where we have used Eq. (O) and 5^ = 2 for the conduc- 
tion band in (100) Si-MOSFET (we have also introduced 
a II to indicate that this is the value of a when the mag- 
netic field is parallel to the 2D electron gas). We assume 
that the only Landau parameter that is crossing an in- 
stability is Fq; the mass ratio m* /nib, contrnlJpd by the 
singlet Landau parameter Ff, is non criticaLclEJ Hence, 
the linear dependence on the density difference n — no of 
the saturation field Hsut measured experimentally follows 
from the Landau phenomenology above. 

Using the experimental data for the i?sat dependence 
on n — no from Refs. g,g and ^ {ni* « 1.5nif,), we find 
that the parameter ajj w 0.6 in the expansion of the Fg 
as a function of (5 = (n — no) /no. 

Although we can explain the linear dependence of -ffsat 
vs. n — no observed for a range of densities in the exper- 
iments, we would like to point out that we cannot easily 
explain the scaling behavior closer to the critical point, as 
reported in Ref. ||. One problem could be that near the 
MIT, with the field vanishing, and with the uncertainty 
in the position of the ferromagnetic transition (notice 
again that in our theory there should be two quantum 
critical points), the exact form for the dependence of the 
saturation field i?sat as a function of T and n may be 



While a parallel magnetic field couples only to the elec- 
trons via the Zeeman coupling, a perpendicular magnetic 
field, H±, also couples to the orbital motion and pro- 
duces Landau levels. In high-density Si-MOSFET sam- 
ples, far from the critical value nc, and for small per- 
pendicular fields, it is observed that the electronic states 
are localized. With increasing magnetic field the elec- 
trons tend to delocafize because of the Lorentz force- 
leading to the so-called fioating of the extended states.cj 
For sufficiently large H±, orbital effects start to play a 
role, and a series qtnuantum Hall phases are observed in 
the Si-MOSFETsa Integer quantum Hall plateaus are 
observed when an integer number of electrons becomes 
commensurate with the number of fiux quanta (Jjq ~ ch/e 
piercing the system. 



H, 



(52) 



where h' is the magnetic filling factor (number of electrons 
per flux quanta). In the case of high-density Si-MOSFET 
the bare Zeeman splitting is very small compared with 
the cyclotron energy. 



huJr 



heHA_ 
nibC 



(53) 



implying that up and down states are nearly degener- 
ate. Moreover, the band structure dispersion contains 
two valleys that are almost degenerate. Equation ( [5^ ) 
for different integers v gives the "Landau fan" in the 
H± X n plane. Whenever condition Eq. (|5l) is obeyed, 
the system sits at the middle of a Hall plateau, where 
the Hall resistivity, p^y, becomes quantized in units of 
h/e'^ {pxy — JiJve^) while the longitudinal resistance, 
Pxx, vanishes.cJ The fact that the quantum Hall effect 
occurs implies that Landau levels are occupied by well- 
defined quasiparticles. This gives extra support to our 
assumption that a Fermi liquid description is correct at 
intermediate temperatures and magnetic fields. 

For a given filling fraction u and at lower densities, 
close to but higher than the zero field nc (see Ref. ^) , the 
quantized Hall states are lost. This occurs at some criti- 
cal density nc{v), at which the states at the center of the 



Landau band localize, and the. longitudinal resistance di- 
verges as T — > 0. It is knownCj that nc{v) varies substan- 
tially along the "Landau fan". In fact, it has been shown 
by Kravchenko and collaborators that if the electron sys- 
tem is spin polarized in the plane (by a parallel mag- 
netic field) and a perpendicular field is applied, rthe lo- 
calization effects occur already at higher densities.ca This 
demonstrates the importance of spin interactions in the 
problem: the correlation effects are reduced when fully 
spin polarization is achieved, making the system closer 
to noninteracting and thus more sensitive to localization 
effects. One of the most interesting observations is that 
for z/ = 4, 8, 12, corresponding to filling factors within the 
cyclotron gaps, the localization occurs at higher densities 
(or magnetic fields fields) than for z/ = 2, 6, 10, where fill- 
ing factors fall into spin gaps (recall that g^ = 2). This 
effect has not been quantitatively explained so far. We 
show below that it can be understood as the localiza- 
tion of the carriers due to the enhancement of the spin 
susceptibility near the zero- field critical density ric. 

As pointed out earlier, at the localization transition 
there is an enhancement of the magnetic susceptibil- 
ity, and consequently of the effective Lande g-factor, g* . 
Therefore, the Zeeman splitting energy 



Ez{n) = g*fiBHA 



(54) 



becomes large and of order of the cyclotron energy. Here 
g* can be obtained from the Fermi liquid theory [see Eq. 
(1)], 



go 



(55) 



which diverges at the quantum critical point. This im- 
plies that there is a level crossing between the down spin 
state of the ith Landau level with the spin up state of the 
i + 1th Landau level. When the crossing occurs there is 
an excess magnetization in the system and, in particular, 
at low field (around 1 Tesla), for the case oi v = A (see 
below), the system becomes fully polarized. The experi- 
mental evidence is that at this point the carriers localize 
since the longitudinal resistivity increases as T ^ O.cil 
This observation implies that below the critical density 
and in the presence of an applied magnetic field, the lo- 
calized state, independent of the direction of the field, is 
indeed spin polarized. 

To illustrate this effect, consider the situation when 
the magnetic field is enough to produce av = A state. In 
this case there are equal number of up and down spins 
filling the first Landau level with energy ?icjc/2. At large 
densities the Zeeman energy is insignificant when com- 
pared to hoJc and the up and down spin states can be 
considered as degenerate. As the density is decreased 
along the curve defined by Eq. (p2) the effective g fac- 
tor increases according to the susceptibility in Eq. (12). 
The first Landau level for up (t) and down (J,) spin states 
changes as a function of density as 



Si=i,t(n) = -hujc - -zEzin) 



(56) 



while the second Landau level changes its energy with 
respect to the density as 



3 1 

Ei=2,^{n) = -huJc - -Ez{n) 

3 1 

Ei=2,i{n) = -huJc + -zEz{n) 



(57) 



Thus, there is a critical density n* such that the previ- 
ously empty second Landau level for spin up, _Ej=2,t('^*)' 
becomes degenerate with the first Landau level with spin 
down, Ei^i^i{n*). Thus, n* is given by 

Ez{,n*)^huo,. (58) 

Using (Q), (p5[) and ^s = eti/{2mQc) we find 



no 1 + 



.90 mt 

2aj_ mo 



(59) 



Using go = 2, nib — 0.2mo, no — 0.8 x 10^^ cm ^, and 
n* = 1.0 X 10" cm-2 (Ref. |4^), we find a± w 0.8. This 
value should be compared with the value of a|| ~ 0.6 
found in the case of a parallel magnetic field. The agree- 
ment is good and gives extra support to the idea that 
the localized state is indeed ferromagnetic even when the 
field is perpendicular to the 2D electron gas. We remark, 
however, that we do not expect the same estimate of a± 
to be applicable to the states i^ = 8 and v = 12, since 
they vanish at densities sufficiently far from Uc to invali- 
date the use of expansion implicit in Eq. (pF 



C. Classical Hall Effect 



A unusual experimental fact related to the behavior of 
the classical Hall coefficient in Si-MOSFET's can also be 
accounted for by our theory. It has been observed that 
the Hall resistance in Si-MOSFET at low temperatures 
and at densities n > ric is insensitive to pajcallel mag- 
netic fields ranging from zero to iJ|[ > i/sat-cH That is, 
there seems to exist a single charge carrier component in 
the metallic phase, for all values of Hsat , instead of two 
independent spin up and down components. This is con- 
sistent with the idea that the conducting fluid present in 
the metallic phase is formed by electron pairs (bosons), 
instead of up/down spin unpaired electrons (fermions). 
Although the number of pairs decreases with an applied 
parallel field (and consequently the number of localized 
electrons increases, see Section HI), notice that the ex- 
periment is performed at fixed current .e3 
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D. A New Experiment: Shot Noise 

The main prediction of our theory is that in the metal- 
Uc phase, at very low temperatures, the system is com- 
posed by incoherent p-wave pairs of electrons. The 
bosons are incoherent because there is no long range or- 
der in the system [see Eq. (|35|)] and therefore no gap in 
the spectrum. This situation is very similar to the case of 
a metallic bosonic liquid.cJ Due to the lack of phase co- 
herence, the usual methods to measure the pair charge, 
such as the Josephson effect, cannot be used. Instead, 
the simplest way to measure the pair charge 2e is by 
making a constriction in the 2D electron density profile 
via external gates and measure the shot noise on the cur- 
rent across the constriction (for shot noise we must have 
ksT ^ eV where V is the voltage applied across the con- 
striction). Since the pairs behave as independent bosonic 
entities, the current fluctuations should be quantized in 
units of the elementary charge 2e.ca 



VI. CONCLUSIONS 

In this paper we argue that one can account, on the 
basis of Landau Fermi liquid theory, for the recent ex- 
perimental observations that the characteristic magnetic 
field Hsat needed to saturate the conductance in 2D 
Si-MOSFETs at low temperatures vanishes at a criti- 
cal value of the electronic density uq. We propose a 
phenomenological expansion for the Landau parameter 
Fq = —l + a{n — no) /no in terms of the electronic den- 
sity n that drives the system ferromagnetic, through a 
Pomeranchuk instability, at ng. As one approaches the 
instability, the spin susceptibility is greatly enhanced, re- 
quiring smaller magnetic fields to fully spin polarize the 
system. At the critical point an arbitrarily small mag- 
netic field fully polarizes the system, since the suscepti- 
bility (at zero temperature) diverges. 

We also analyze the effects of a perpendicular magnetic 
field through the system, in the quantized Hall regime, 
and show that the critical density for the localization of 
the u = A state as compared to the u = 2 state can be ac- 
counted for by considering the crossover between a spin 
polarized and unpolarized state due to the enhancement 
of the Lande g factor. The values of the parameter a es- 
timated separately from parallel and perpendicular field 
experiments agree within 20%. 

In the paramagnetic side, but close to the instability, 
the enhanced spin fluctuations can lead to an attractive 
interaction in the spin triplet channel, similarly to su- 
perfluid "^He. We analyze a Landau-Ginzburg mean-field 
theory that combines p-wave superconductivity and fer- 
romagnetism, and find two quantum critical points as a 
function of density, tiq where ferromagnetism begins, and 
nj, where p-wave pairing ceases. There is an intermediate 
range of densities where p-wave pairing and ferromag- 
netism coexist. In this range, the p-wave state is in a 



nonunitary phase. 

For large enough densities, above a value nA, the para- 
magnon exchange mechanism responsible for attractive 
interactions should cease (when, for example, Fq > 0). 
At these densities, an Anderson insulating state, simi- 
larly to the case of noninteracting electrons, should oc- 
cur. However, the presence of the p-wave paired state for 
Uc < n < nA does not rule out the possibility of a con- 
ducting phase. Since the order parameter for the p-wave 
state is a vector, no order (even algebraic) exists at finite 
temperature, and true superconductivity should only oc- 
cur at T = 0. In this paper we do not present any ex- 
planation why the p-wave pairing would lead to the con- 
ducting phase at finite T; however, we strongly believe 
that if this correlated state exists in the 2D electronic 
systems, it may provide the origin of the extended state 
that continues to accumulate experimental support .e2I 

Summarizing, we propose that the metallic state close 
to the metal-insulator transition in the 2D electron gas 
problem is due to the existence of a paired p-wave state 
close to a ferromagnetic insulating phase. The pairing is 
generated by long wavelength magnetic fluctuations close 
to the quantum critical point. We describe the pairing 
within Landau's Fermi liquid phenomenology and show 
that it provides a consistent description of the data for 
parallel and perpendicular magnetic fields. Moreover, we 
propose new shot noise experiments that can test our 
theory. 
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